Controlling the spontaneous spiking regularity via channel 
blocking on Newman- Watts networks of Hodgkin-Huxley neurons 
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Abstract. - We investigate the regularity of spontaneous spiking activity on Newman- Watts 
small-world networks consisting of biophysically realistic Hodgkin-Huxley neurons with a tunable 
intensity of intrinsic noise and fraction of blocked voltage-gated sodium and potassium ion channels 
embedded in neuronal membranes. We show that there exists an optimal fraction of shortcut 
links between physically distant neurons, as well as an optimal intensity of intrinsic noise, which 
warrant an optimally ordered spontaneous spiking activity. This doubly coherence resonance-like 
phenomenon depends significantly, and can be controlled via the fraction of closed sodium and 
potassium ion channels, whereby the impacts can be understood via the analysis of the firing rate 
function as well as the deterministic system dynamics. Potential biological implications of our 
findings for information propagation across neural networks are also discussed. 



Introduction. — Information processing within the 
nervous system involves noisy components [1]. One ma- 
jor source of noise within neurons is due to the stochas- 
tic dynamics of voltage-gated ion channels embedded in 
neuronal membranes [2]. Voltage-gated ion channels are 
involved in generating and propagating electrical signals 
through neuronal membranes. Hodgkin and Huxley (HH) 
[3] first introduced conductance-based models of these 
channels. However, the original HH model does not 
take into account the stochastic dynamics of voltage-gated 
channels, but provides a deterministic description of the 
membrane potential, which is valid only within the limit 
of very large cell sizes. When the population of ion chan- 
nels is finite, the stochastic dynamics of voltage-gated ion 
channels (or ion channel noise) causes subthreshold fluc- 
tuations in the membrane voltage [2] , and can have signif- 
icant impacts on the neuronal dynamics [4-11]. 

The intensity of the channel noise is related to the num- 
ber of ion channels, but its actual impact is determined by 
the number of channels that are open near the threshold 
for spike firing. Schncidman et al. [6] showed that there 
is a short distance in terms of the number of open chan- 
nels between firing and non-firing stable states, and that 



fluctuations due to only a few channels are responsible 
for the transition between these two stable states. They 
also reported that this spontaneous transition is the cause 
for the missing spikes, the subthreshold membrane volt- 
age, as well as spontaneous spikes. Therefore, controlling 
the number of working ion channels for a given membrane 
patch is of great importance, particularly to understand 
the impact of a specific ion channel type on neuronal dy- 
namics. In this context, some toxins such as tetraethylam- 
monium (TEA), tetradotoxin (TTX) and saxitoxin (STX) 
are widely used in experiments to block or reduce the 
number of specific ion channels [12]. It is also possible 
to examine the effects of changing the number of specific 
ion channels on neuronal dynamics through computational 
models. Schmid et al. [13,14] investigated the regularity of 
spontaneous spiking activity of a stochastic HH model at a 
single neuron level, and showed that it is possible to either 
increase or decrease the regularity of spontaneous spike 
trains by blocking some portion of either potassium or 
sodium ion channels. Recently, Gong et al. [15] extended 
the subject by examining the effects of channel blockage 
on the collective spontaneous spiking activity of coupled 
stochastic HH neurons, and found that the latter can be 
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reduced or enhanced by blocking the sodium or potassium 
channels within small or large cell sizes. They constructed 
the interaction network as an array of bi-directionally cou- 
pled neurons. However, small-world (SW) networks, com- 
bining high clustering with scattered long range connec- 
tions [16] have been suggested as an attractive option 
to search for connectivity information of both anatomi- 
cal and functional networks in the brain, mainly because 
this topology can support both local and distributed in- 
formation processing [17]. Therefore, SW networks have 
been widely used to understand how neuronal circuitry 
generates complex patterns of activity [18-22]. 

Our aim in this Letter is to extend the subject by us- 
ing SW networks as the underlying interaction topology 
between neurons, and to investigate the impact of sodium 
and potassium channel blockage on the spontaneous col- 
lective spiking regularity as a function of the network 
structure and cell size in a manner that more closely mim- 
ics actual conditions. 

Mathematical model and setup. Within this 
study we use the HH model of neuronal dynamics as the 
basic building block of the investigated system, whereby 
the time evolution of the membrane potential for coupled 
neurons on an arbitrary network in the absence of deter- 
ministic external stimuli is given as follows: 

C m dV t /dt = G Na (m t , hi){V Na - Vi) + G K (n t ){V K - V % ) 
+G L (V L -V i )+Y,Zij(V j -V i ) (1) 

3 

Here Vi denotes the membrane potential of neuron i = 
1,...,N (N being the system size), C m = l^Fcm -2 is 
the membrane capacity, and Gjv a , Gk and Gl repre- 
sent sodium, potassium and leakage conductances, respec- 
tively. V Na = 50mV, V K = -77mV and V L = -54.4mV 
are the reversal potentials for the sodium, potassium and 
leakage channels, respectively. Moreover, denotes the 
coupling strength between neurons i and j , whereby we set 
Eij = s if the two are connected or £y = otherwise. In 
the model, the leakage conductance is assumed to be con- 
stant, equaling Gl = 0.3mScm~ 2 , while the sodium and 
potassium conductances change dynamically according to 
the following two equations: 

G Na {m u hi) = g^XNamfhi, G K {m) = tfF*x K n\ (2) 

Here g^ x = 120mSc m - 2 and g% ax = 36mSc m - 2 are 
the maximal sodium and potassium conductances, respec- 
tively. Moreover, m and h denote the activation and inac- 
tivation gating variables for the sodium channel, respec- 
tively, whereas the potassium channel includes an activa- 
tion gating variable n. In Eq. (2) we also introduce two 
scaling factors, xjy a and xk, which are the fractions of 
non-blocked ion channels to the total number of sodium 
{Nno,) or potassium (Nk) ion channels within the patch 
area, respectively [13-15]. These factors are confined to 
the unit interval. 



In the HH model, activation and inactivation gating 
variables, rrii, i%i and hi, change over time in response 
to the membrane potential following first-order differential 
equations within the limit of very large cell sizes. However, 
when the population of ion channels is finite, the stochas- 
tic behavior of voltage-gated ion channels (or ion channel 
noise) must be taken appropriately into account. In this 
study, we follow the approach in previous works related to 
the channel block [13-15], and use the algorithm proposed 
by Fox [23]. In the Fox's algorithm, variables of stochas- 
tic gating dynamics are described with the corresponding 
Langevin generalization 

dxi/dt = a x (l — Xi) — P x Xi + £ Xi (t), Xi=mi,m,hi (3) 

where a x and (3 X are rate functions for the gating variable 
Xi. The probabilistic nature of the channels appears as a 
source of noise £, Xi (t) in Eq. (3), which is an independent 
zero mean Gaussian noise whose autocorrelation function 
is given as follows [23]: 



<fm(*)£m(f)> = 



<&(*)&(*')> = 



2c^ m /3 r7 



N Na x Na (a m + (3 m ) 



S(t-t') (4) 



N Na x Na (a h + Ph) 



6(t-t>) (5) 



2a„/3„ 



N K x K {a n + Pr, 



-S(t-t') 



(6) 



where the factors, xj\r a and xk , are used again to disre- 
gard the blocked channels, which do not contribute to the 
intrinsic channel noise. Given the assumption of homoge- 
neous sodium and potassium ion channel densities, chan- 
nel numbers are calculated via iVjv a = PNaS, Nk = pkS 
where p^ a = 60^m -2 and px = 18/im -2 are the sodium 
and potassium channel densities, respectively, whereas S 
represents the total membrane area or the cell size [23]. 
Equations (l)-(6) constitute the stochastic HH model, 
where the cell size S determines the intrinsic noise level 
through the number of ion channels N^ a and Nk- When 
S is large stochastic effects related to the channel noise 
are negligible due to large numbers of ion channels in Eqs. 
(4)- (6), and thus the intrinsic channel noise appearing in 
Eq. (3) for the gating variables vanishes. Accordingly, 
the stochastic model then approaches the deterministic 
description. However, when the number of ion channels 
(or the cell size S) is small, stochastic effects have a sig- 
nificant impact on the membrane dynamics [13-15,23,24]. 

Furthermore, the interaction network is comprised of 
identical HH neurons, initially each having connectivity 
k = 2, with the system size set to N — 60. Following the 
Newman- Watts model [25], we start with a regular ring 
[see Fig. 1(a)] and make a random draw of two neurons. 
Subsequently, if they are not already connected, we add a 
non-directed link between them. This process is repeated 
until a total of M new links are added [see Fig. 1(b)], 
finally resulting in a network which is equivalent to the 
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Fig. 1: Examples of considered network topologies. For clarity 
regarding k and p only N = 25 vertices are displayed in each 
panel, (a) Regular ring characterized by p — with periodic 
boundary conditions. Each vertex is connected to its k = 2 
nearest neighbors, (b) Realization of Newman-Watts small- 
world topology via randomly adding M new links (in this case 
6 new links were added, hence p — 0.02). 




one obtained if new edges would be added with probability 
p = 2M/[N(N- 1)]. 

In order to asses the system's dynamics, we measured 
the collective temporal behavior of the network by first 
calculating the average membrane potential, V avg {t) = 
N^ 1 ^ ji=1 Vi{t) } corresponding to the mean field of the 
network, and then quantifying the spontaneous collective 
spiking regularity by the coefficient of variation (CV) of 
the inter-spike intervals (I Sis) according to: 

cv -JTsT)- (isT) (7) 

Here 'spike' times were defined by the upward crossing 
of V avg (t) past a detection threshold of OmV (referred to 
below as 'network spikes'), whereas (ISI) and (ISI 2 ) de- 
note the mean and the mean squared inter-spike intervals, 
respectively. Notably, CV is characterized by smaller val- 
ues for more ordered spike trains and vanishes for a de- 
terministic signal. In the results presented below, for any 
given set of the remaining network parameters (scaling fac- 
tor, noise level, coupling strength), quantitative results re 
pooled from simulations of 10 realizations of the network 
for any given value of p. 

Results and discussion. — In what follows, we sys- 
temically analyze the impact of sodium and potassium ion 
channel block on the spontaneous collective spiking regu- 
larity of the network as a function of the network topology 
and cell size (thus, ion channel noise). First, we investi- 
gate how CV changes with the scaling factors, xjsr a or xk, 
as a function of the network topology by a fixed cell size 
S = 6/im 2 and coupling strength e = 0.1. Thereby we 
vary the density of one channel type (either x^ a or xk) 
while keeping the other equal to one. Obtained results 
are presented in Fig. 2(a) for sodium and in Fig. 2(b) for 
potassium ion channel block. Notably, we did not consider 
XNa < 0.8 because then the average membrane potential 
V aV g did not include spikes. For all values of XNa and 
for all xk > 0.5 considered in Fig. 2, the CV exhibits a 
well-expressed minimum by an optimal fraction of random 
shortcuts p, at which the spontaneous collective spiking is 
most regular. Moreover, it is evident that there exists an 
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Fig. 2: The dependence of the spontaneous spiking regularity 
(CV) on p obtained by a fixed patch area S = 6/xm 2 and cou- 
pling strength e = 0.1. (a) Different levels of sodium channel 
block, (b) Different levels of potassium channel block. 

optimal p independent of the scaling factors, thus indicat- 
ing the robust existence of topology-dependent coherence 
resonance [26-29]. Interestingly, an optimal topology is 
obtained at p ~ 0.15 for both the sodium as well as the 
potassium channel blocks, indicating the existence of a 
universally optimal network structure. However, as xk 
decreases to xk < 0.5, the minimal values of CV arc ob- 
tained already by slightly lower values of p, yet after reach- 
ing the minimum, the regularity of V avg becomes constant 
and virtually independent of the network topology. In- 
deed, the level of Xk in the range < 0.5 plays a rather in- 
significant role for the temporal regularity of the average 
membrane potential for all p > 0.075, as can be inferred 
from the bottom two curves depicted in Fig. 2(b). 

Furthermore, it is worth noting that the reduction of 
working sodium channels [decrease of XNa in Fig. 2(a)] 
decreases the collective spiking regularity, whereas a re- 
duction of working potassium channels [decrease of xk in 
Fig. 2(b)] increases it by a fixed cell size, as reported also 
in [13, 14] for a single neuron. Notably, Gong et al. [15] 
found that sodium or potassium channel blocks can either 
enhance or reduce the collective spiking regularity of an ar- 
ray of bi-directionally coupled neurons, depending on the 
cell size. In our case, however, for p w 0.15 the collective 
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Fig. 3: The dependence of the mean firing rate (Fr) on p 
obtained by a fixed patch area S = and coupling strength 

e = 0.1. (a) Different levels of sodium channel block, (b) 
Different levels of potassium channel block. 

regularity is very high for all values of xk and xjq a > 0.9 
if compared to the values reported in [15] for the same cell 
size and scaling factors (see Figs. 1 and 2 in [15]). This 
indeed strongly supports the fact that the fine-tuning of 
the network structure via p is able to significantly facili- 
tate the spontaneous spiking regularity of neurons in the 
presence of either sodium or potassium ion channel block. 

To gain more insight into the dependence of the col- 
lective spiking regularity on p and the scaling factors, we 
further calculate the mean firing rate (Fr) for different val- 
ues of p and XNa or xk- Obtained results are presented 
in Fig. 3(a) for sodium and in Fig. 3(b) for potassium ion 
channel block. In case of sodium channel block, the mean 
firing rate exhibits a resonance-like behavior with a dis- 
tinct maximum appearing by the optimal fraction of ran- 
dom shortcuts equaling approximately p w 0.15. Potas- 
sium channel block leads to similar behavior, but only 
when xk — ► 1. As xk decreases below 0.95, the maxima 
are achieved by fewer added shortcuts (lower values of p). 
Thus, the outlay of CV presented in Fig. 2 follows the 
characteristics of the firing rate, and moreover, decreases 
with higher Fr (compare Figs. 2 and 3). In this sense, 
results presented in Fig. 3 are consistent with those pre- 
sented in Fig. 2 for both types of ion channel blockage. In 



particular, the refractory period of spike generation postu- 
lates a certain minimal time that needs to elapse between 
two spikes. Thus, as the firing frequency increases, the 
mean ISI decreases and can approach the refractory pe- 
riod [14]. In other words, as the firing frequency increases, 
the refractory period becomes more significant relative to 
the diminishing ISI, contributing to the reduction of the 
output CV [30,31]. We therefore argue that the network 
topology affects the regularity of spontaneous neuronal 
spiking by means of serving as a scaling factor for the 
firing rate. 

Related to the above argumentation, it is interesting to 
elaborate further on the results obtained for xk < 0.5. 
There a further increase in p, past the optimal value for 
which the maximal firing rate is obtained, decreases the 
mean firing rate fairly slightly [see the upper two curves 
in Fig. 3(b)]. As shown in the bifurcation diagram for 
the deterministic autonomous HH model in [14] [see their 
Fig. 1(b)], with decreasing the potassium conductance 
via xk in Eq. (2), a sub-critical Hopf bifurcation occurs 
at xk = 0.549, after which the system enters a region of 
stable oscillatory spiking solutions within the purely de- 
terministic description [13,14]. We calculated the firing 
rate of the deterministic network by considering just the 
conductance reductions xk = 0.5 and xk = 0.25, and 
found them equaling 51Hz and 61Hz, respectively. These 
values are almost equal to those obtained for the model 
with stochastic components as p approaches to its upper 
boundary (all-to-all connected network) [see Fig. 3(b)]. 
This suggests that the stochastic dynamics has little effect 
on the spontaneous regularity for the potassium channel 
block if xk < 0.5. It also suggests that the network con- 
nectivity has a negligible impact on the firing rate (or CV) 
for all xk < 0.5 as p — » 1. Importantly however, for small 
values of p the firing rate reaches a maximum above its 
deterministic value, indicating the existence of an optimal 
small-world topology that warrants the best regularity of 
spontaneous collective spiking activity even in a region 
with a purely deterministic description, which then again 
disappears as p — > (regular nearest-neighbor network). 

Following above investigations, where we found that 
the CV exhibits a well-expressed minimum at an optimal 
p « 0.15 by a fixed cell size S = 6/zm 2 , we investigate now 
how the regularity by the optimal value of p changes with 
the cell size for different values of the scaling factors de- 
termining the non-blocked channel densities x^a and xk- 
Obtained results are presented in Fig. 4(a) for sodium 
and in Fig. 4(b) for potassium ion channel block. We find 
that the spontaneous collective spiking regularity depends 
strongly on the cell size for all values of XMa and for all 
xk > 0.5, with the best regularity emerging for small, but 
non-zero cell sizes of around S = 2 — 4/im 2 constituting 
a specific intensity of intrinsic noise. For this range of S 
values, the regularity changes slightly depending on xk, 
whereas it may change rather substantially depending on 
XNa- However, for larger cells (large S) the CV differs 
slightly depending on XNa while it differs substantially 
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Fig. 4: The dependence of the spontaneous spiking regularity 
(CV) on the patch area S obtained by the coupling strength 
e = 0.1 and p — 0.15. (a) Different levels of sodium channel 
block, (b) Different levels of potassium channel block. 

depending on xk- 

On the other hand, for xk < 0.5 an increase in the 
cell size beyond a value of around S = 2p,m 2 almost does 
not improve the regularity of the output, albeit a small 
improvement in CV can be inferred for xk = 0.25 if com- 
pared to xk = 0.5. This collective behavior seems similar 
to that of a single HH neuron reported in [13, 14], where 
the spontaneous regularity was found to increase as the 
cell size increased for xk < 0.5. However, we show here 
that this increase is larger for a single neuron if compared 
to that of the whole network. Also, we note that the 
behavior of an array of bi-directionally coupled neurons 
reported in [15] was found to be rather different, where 
indeed the regularity decreased as the cell size increased 
past S > 2pm 2 . 

Notably, we also investigated the CV dependencies pre- 
sented in Fig. 4 via the mean firing rate (not shown), 
and arrived at the findings that are in agreement with 
the above-established link between the ISI and the fir- 
ing rate when interpreting results presented in Figs. 2 
and 3. The latter support the fact that, not only does 
the network topology affect the regularity of spontaneous 
neuronal spiking by means of serving as a scaling factor 
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Fig. 5: The dependence of the spontaneous spiking regularity 
(CV) on p for three different coupling strengths e obtained by 
a fixed patch area S = 6/xm 2 . (a) Sodium channel block is 
XNa — 0.95. (b) Potassium channel block is xk — 0.85. 

for the firing rate, but the regularity by the optimal SW 
topology is affected also by the intensity of intrinsic noise, 
similarly as reported earlier for the internal-noise coher- 
ence resonance. These two facts combined yield a doubly 
coherence resonance-like phenomenon [32,33] in our case, 
related to the regularity of spontaneous neuronal spiking, 
that is brought about by the fine-tuning of both p and S. 

In order to further extend and conclude the study, we 
investigate the dependence of the spontaneous collective 
regularity on p for three different coupling strengths e, by 
a fixed cell size S = 6^m 2 , and two different cases of ion 
channel block x^ a = 0.95 and xk = 0.85. The two values 
were selected so as to lead to a coherence resonance in CV 
against p. Obtained results are presented in Fig. 5. From 
them we arrive at two findings that are independent of the 
type of the ion channel block. First, the resonance- like 
outlay of CV in dependence on p is similar for all cou- 
pling strengths and its overall minimum does not change. 
Second, the coupling strength is inversely correlated with 
the value of p that yields the minimum CV. Accordingly, 
as e increases the minimum CV is obtained by fewer and 
fewer added shortcuts (by lower values of p). However, 
these correlations are not completely valid for the total 
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span of p, because irrespective of e, the spontaneous col- 
lective regularity becomes very low for extreme values of p 
(i.e. p — ► or p — » 1). Nevertheless, the findings are con- 
sistent with previous results related to different subjects 
but to the same network structure [22,34]. 

Summary. — In sum, we show that the regularity of 
spontaneous spiking activity can be resonantly enhanced 
via fine-tuning of both p and S. The extend of improved 
regularity is thereby maximal by high XNa and low xk 
values. Solely by potassium ion channel blockage consti- 
tuted by xk < 0.5 the resonant dependence on p van- 
ishes and the first occurrence of the minimal CV occurs 
already by somewhat lower values of p. There also ex- 
ists an optimal intensity of internal noise, constituted by 
cell sizes spanning S = 2 — 4/im 2 , by which the regularity 
of spontaneous spiking activity on small-world networks 
is maximal. Such cell sizes may correspond to so-called 
'hot spots' [35] or otherwise spatially restricted regions in 
dendritic or axonal trees [36,37]. Both resonant-like de- 
pendencies by different values of x N a and Xk can be nicely 
corroborated by the analysis of the firing rate as well as 
the deterministic model dynamics. Finally, we show that 
the resonant outlay of the CV dependence on p prevails 
irrespective of the coupling strength, although the minima 
shift in accordance with some previous findings. Although 
the present model considers the intrinsic noise stemming 
from the stochastic dynamics of ion channels, our study 
also provides some indication towards the potential im- 
pact of synaptic noise because the latter may be approxi- 
mated by a non-zero mean Gaussian white noise [38] . We 
hope our study will prove useful in striving towards the 
understanding of importance of structure [39-42] and un- 
certainty in active neuronal networks. 
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